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Abstract 

If a body of dielectric material is coated by a plasmonic structure of negative dielectric 
material with nonzero loss parameter, then cloaking by anomalous localized resonance (CALR) 
may occur as the loss parameter tends to zero. It was proved in [1] [2] that if the coated 
structure is circular (2D) and dielectric constant of the shell is a negative constant (with loss 
parameter), then CALR occurs, and if the coated structure is spherical (3D), then CALR does 
not occur. The aim of this paper is to show that the CALR takes place if the spherical coated 
structure has a specially designed anisotropic dielectric tensor. The anisotropic dielectric 
tensor is designed by unfolding a folded geometry. 

1 Introduction 

If a body of dielectric material (core) is coated by a plasmonic structure of negative dielectric 
constant with nonzero loss parameter (shell), then anomalous localized resonance may occur as 
the loss parameter tends to zero. To be precise, let f2 be a bounded domain in Mf^, d = 2,3, and D 
be a domain whose closure is contained in $7. In other words, D is the core and \ D is the shell. 
For a given loss parameter S > 0, the permittivity distribution in is given by 

in M'^ \ n, 

es={es + i6 mn\D, (1.1) 
in D. 

Here ec is a positive constant, but is a negative constant representing the negative dielectric 
constant of the shell. For a given function / compactly supported in \ fl satisfying 

fdx = (1.2) 

(which is required by conservation of charge), we consider the following dielectric problem: 

V-esVVs = f iuR'*, (1.3) 
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with the decay condition Vs{x) — >■ as |a;| — > cxj. The equation (|1.3[) is known as the quasistatic 
equation and the real part of — VV^ (x)e~"^*, where w is the frequency and t is the time, represents 
an approximation for the physical electric field in the vicinity of il, when the wavelength of the 
electromagnetic radiation is large compared to f2. 
Let 

S,- :=3 / ealVFif dx= / 5\VVs\^ dx (1.4) 

Js.'i jn\D 

(3 for the imaginary part), which, within a factor proportional to the frequency, approximately 
represents the time averaged electromagnetic power produced by the source dissipated into heat. 
Also for any region T let 

i?°(T)= / \\7Vs\'dx (1.5) 

JT 

which when T is outside fl approximately represents, within a proportionality constant, the time 
averaged electrical energy stored in the region T. Anomalous localized resonance is the phenomenon 
of field blow-up in a localized region. It may (and may not) occur depending upon the structure 
and the location of the source. Quantitatively, it is characterized by Eg{T) — > oo as (5 — ?> for all 
regions T that overlap the region of anomalous resonance, and this defines that region. Cloaking 
due to anomalous localized resonance (CALR) may occur when the support of the source, or part 
of it, lies in the anomalously resonant region. Physically the enormous fields in the anomalously 
resonant region interact with the source to create a sort of optical molasses, against which the 
source has to do a tremendous amount of work to maintain its amplitude, and this work tends to 
infinity as (5 — )■ 0. Quantitatively it is characterized by iJ^ — > oo as 5 — >■ oo. 

This phenomena of anomolous resonance was first discovered by Nicorovici, McPhcdran and 
Milton [TS] and is related to invisibility cloaking [TT]: the localized resonant fields created by a 
source can act back on the source and mask it (assuming the source is normalized to produce fixed 
power). It is also related to superlenses [TBI [17] since, as shown in [15], the anomalous resonance 
can create apparent point sources. For these connections and further developments tied to this 
form of invisibility cloaking, we refer to [U [U |31 HI [IH] and references therein. Anomalous resonance 
is also presumably responsible for cloaking due to complementary media (HI [IHl [H] , although we 
do not study this here. 

The problem of cloaking by anomalous localized resonance (CALR) can be formulated as the 
problem of identifying the sources / such that first 

Es:=l S\S/Vs\^ dx ^ oo as S ^ 0, (1.6) 
Jn\D 

and secondly, Vs/\/Es goes to zero outside some radius a, as (5 — > 0: 

|V5(a;)/V^I ^ as 5^0 when|a;|>a. (1.7) 

Since the quantity Eg is proportional to the electromagnetic power dissipated into heat by the time 
harmonic electrical field averaged over time, (|1.6p implies an infinite amount of energy dissipated 
per unit time in the limit (5 — >■ which is unphysical. If we rescale the source / by a factor of 
1/\/Es then the source will produce the same power independently of S and the new associated 
potential Vs/^/Es will, by (jl.7[) . approach zero outside the radius a. Hence, cloaking due to 
anomalous localized resonance (CALR) occurs. The normalized source is essentially invisible from 
the outside, yet the fields inside are very large. 

In the recent papers [TJ [2] the authors developed a spectral approach to analyze the CALR 
phenomenon. In particular, they show that if D and are concentric disks in of radii and r^, 
respectively, and ~ —1, then there is a critical radius r» such that for any source / supported 
outside CALR docs not occur, and for sources / satisfying a mild (gap) condition CALR takes 
place. The critical radius is given by = i/rf /r^ if = 1, and by r^, = rl/vi if ec 7^ 1. It is also 
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proved that if 7^ —1, then CALR does not occur: Es is bounded regardless of S and the location 
of the source. It is worth mentioning that these results (when Ec = —eg = 1) were extended in [7] 
to the case when the core D is not radial by a different method based on a variational approach. 
There the source / is assumed to be supported on circles. 

The situation in three dimensions is completely different. If D and fl arc concentric balls in M"^, 
CALR does not occur whatever and ec are, as long as they are constants. We emphasize that 
this discrepancy comes from the convergence rate of the singular values of the Neumann-Poincare- 
type operator associated with the structure. In 2D, they converge to exponentially fast, but in 
3D they converge only at the rate of 1/n. See [2]. The absence of CALR in such coated sphere 
geometries is also linked with the absence of perfect plasmon waves: see the appendix in [7]. On 
the other hand, in a slab geometry CALR is known to occur in three dimensions with a single 
dipolar source [11]. (CALR is also known to occur for the full time-harmonic Maxwell equations 
with a single dipolar source outside the slab superlens [51 [TTl [T5] .) 

The purpose of this paper is to show that we are able to make CALR occur in three dimensions 
by using a shell with a specially designed anisotropic dielectric constant. In fact, let D and fl be 
concentric balls in R"^ of radii and , and choose tq so that ro > re . For a given loss parameter 
S > 0, define the dielectric constant by 



e5(x) 



6(5-2|xn 




|x| > r-e, 

ri < |x| < re, 

|x| < ri, 



(1.8) 



where I is the 3x3 identity matrix, and Ec constants, x 



and 



>0, 



(1 + a)re 



(1.9) 



Note that is anisotropic and variable in the shell. This dielectric constant is obtained by push- 
forwarding (unfolding) that of a folded geometry as in Figured] (See the next section for details.) It 
is worth mentioning that this idea of a folded geometry has been used in [12] to prove CALR in the 
analogous two-dimensional cylinder structure for a finite set of dipolar sources. Folded geometries 
were first introduced in [S] to explain the properties of super lenses, and their unfolding map was 
generalized in [12] to allow for three different fields, rather than a single one, in the overlapping 
regions. Folded cylinder structures were studied as superlenses in [20] and folded geometries using 
bipolar coordinates were introduced in to obtain new complementary media cloaking structures. 
More general folded geometries were rigorously investigated in |14j . 



* 





















re ro 



Figure 1: unfolding map 
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For a given source / supported outside Br^ let Vs be the solution to 

(\/ ■ {esVVs) ^ f iuR^ ^^^^^^ 
\Vs{x) ^ as |x| oo, 

and define 

Es^'i [ esVVs-VVsd^, (1.11) 
where Vg is the complex conjugate of Vs. Let F be the Newtonian potential of the source /, i.e., 

F(x) / G(x-y)/(y)dy, (1.12) 

with G'(x — y) = — 4^|x-y| ■ Since / is supported in M.^ \ Br^, F is harmonic in |x| < R for some 
R> re and can be expressed there as 

oo n 

^W=E E /nN"^n(x), (1.13) 

Ti=0 k= — n 

where ^^(x) is the (real) spherical harmonic of degree n and order k. 
The following is the main result of this paper. 

Theorem 1.1 (i) If Ec ~ ~ I7 then weak CALR occurs and the critical radius is r* = %Jr^rQ, 
i.e., if the source function f is supported inside the sphere of radius r* (and its Newtonian 
potential does not extend harmonically to R^J, then 



limsup£'5 = 00, (1-14) 

<5-!-0 

and there exists a constant C such that 

\V5{^)\ < C (1.15) 

for all X with |x| > r^r^^^ . If, in addition, the Fourier coefficients of F satisfy the 
following gap condition: 

[GCl].' There exists a sequence {rij} with Tii < 712 < • • • such that 

J— >CXD ^ ^ 

k——nj 

where p :— r^/rQ, then CALR occurs, i.e., 

lim Es = 00, (1.16) 

<5-i-0 

and Vs/\/Es goes to zero outside the radius r^/r^, as implied by U.IS]) . 

(a) If ec 7^ — = 1, then weak CALR occurs and the critical radius is — rg. //, in addition, 
the Fourier coefficients f^^ of F satisfy 

[GC2]; There exists a sequence {uj} with ni < n2 < ■ ■ ■ such that 

k— — nj 
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then CALR occurs. 

(Hi) If ~es 7^ 1; then CALR does not occur. 

We emphasize that [GCl] and [GC2] are mild conditions on the Fourier coefficients of the 
Newtonian potential of the source function. For example, if the source function is a dipole in 
Br, \Be, i.e., /(x) = a- V(5y(x) for a vector a and y £ Br, \i?e where Sy is the Dirac delta function 
at y. [GCl] and [GC2] hold and CALR takes place. A proof of this fact is provided in the appendix. 
Similarly one can show that if / is a quadrupole, f{x) = A : VV(5y(x) = ^^-^ Uij q^q^ 5y (x) for 
a 3 X 3 matrix A = (a,j) and y G Br, \ B^, then [GCl] and [GC2] hold. 



2 Proof of Theorem 11.11 



Let ri, Te and 7'o be positive constants satisfying r; < re < rg, as before. For given constants Kc 
Kg and Km, and a source / supported outside Br^, let Uc, Ug and Um be the functions satisfying 

Auc = in Br„ , 



Aus = in Bro \ Br^ , 
Au„t = .f in R3 \ R;r^ 





duc 


dus 


Uc = Us, 


or 


or 




dus 


dUm 


Us — '^7n-) 


Or 


or 


u„(x) 


as [x 


— ?> CO. 



on dBrg, 
on dBr 



(2.1) 



We emphasize that the domains of Uc, Ug, and Um are overlapping on < [x[ < ro so that the 
solutions combined may be considered as the solution of the transmission problem with dielectric 
constants Kc, and k„i in the folded geometry as shown in Figure [1] We unfold the solution into 
one whose domain is not overlapping, following the idea in |12] . 

In terms of spherical coordinates (r, 9, </>), the unfolding map $ = {^c, ^s, ^m} is given by 



f$„(r,0,0) = (r,0,0). 



ra - re 



{r - re),( 



-r,e,^ 

Jo 



r > re, 

re <r < ro, 

r < ro. 



(2.2) 



Then the folding map can be written (with an abuse of notation) as 

|x| > re, 
rj < [x[ < re, 
|x| < r„ 

where a and b are constants defined in (jl.9|) . 

Let /«(x) be the push-forward by the unfolding map $, namely. 




(2.3) 



/t(x) 



'K„[detV$™(y)[-iV$„,(y)V$„(y)*, [x[ > re, 
-«;,[detV$,(y)[-iV$,(y)V$,(y)*, r, < [x[ < re, 

Ke[detV$e(y)|-^V$e(y)V$e(y)*, |x[ < r„ 



(2.4) 



where x = *&(y)- By straight-forward computations one can see that n = e given in (|1.8[) if we set 

Kjn = 1, Kg = -{eg + iS), Ke = Cc- (2.5) 
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Moreover, the solution Vs to (jl.lO[) is given by 



\t„o$-i(x) if|x|>re, 
\/5(x) = <( o $-i(x) ifrj<|x|<re, 
Uc o $~i(x) if |x| < ri, 



(2.6) 



and by the change of variables x = ^sCy) and (|2.4p . we have 
Es^'i [ e{x)\/Vs{x)-VVs{x)^S [ 

JK3 Jr, 



<|y|<'"o 



|V^^,s(y)r 



(2.7) 



Suppose that the source / is supported in |x| > R for some R > Te- Then the solution u to 
(|2.ip can be expressed in |x| < i? as follows: 



oo n 



^^W = E E if|x|<ro, 

n— k——n 
oo n 

= E E (^"1^1" + c'M-"-')Yj:{^), if r, < |x| < ro, 

n—0 k——n 
oo n 

i„(x) -EE (^ni^i" + rff;ixr"-^)i;f (X), if ^ < ix| < r, 

n—O k——n 

where the coefficients satisfy the foUowing relations resulting from the interface conditions: 



(2.8) 

















n-1 








= K.(6!;nrJ - 


cf,(n + 




), 




= '*m(eJ^?T.rg - 


df,(n 







By solving this system of linear equations one can see that 

„k _ k uk _ 7 k k _ k jk _ , k 

where 

„2n+l 



p2"+l(2n + l)2K^K, 



„2ri+l 



-p^"+V,„(2n + l)((n + 1)ks + nUc) 



(77.2 _|_ yj)(K3 - Kc)(Ks - «:m) " p2"+H('^ + + nKc){{n + l)Kjn + TlKs) ' 

^2n+l 



{■n? + n)(Ks - Kc){ks - Km) - /92"+i((n + 1)ks + nKc){{n + 1)k„i + uKs) ' 
^ _ nrg"+^[p2"+^(K„ - ^^)(("- + 1)ks + »^c) + (ks - Kc){nKm + (n + 1)k^)] 

(n2 + n){K.s - Kc)(Ks - Km) - p2n+l((7j + 1)^^ + nKc){{n + 1)Kto + UKs) ' 

Here p is defined to be Te/ro 

Let F be the Newtonian potential of /, as before. Since m — F is harmonic in |x| 
tends to as |x| — > oo, we have 

k rk 
'^n Jn- 

So Um (the solution in the matrix) is given by 



Umi^) = ^^(X) + E E /nrf«|x| — (X). 



(2.9) 
(2.10) 
(2.11) 
(2.12) 

> and 
(2.13) 

(2.14) 



n=0 k=—n 
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Since \dn\ < Cr^ , we have 

oo n 

|u„(x) -F(x)| <CY,Y. l/nko"|x|— ^ < oo (2.15) 



Ti=0 k=~n 



if |x| = r > rpTe ^. This proves (jl.l5|) . 
The solution in the shell is given by 



oo n 



"^(y) = E E /„'(^>«iyr+c„iyr"-^)r,f(y). (2.16) 

n— k— — n 

Using Green's identity and the orthogonality of Y^^, we obtain that 



dus f dug 
or 



00 n 

n=0 k— — n 

00 n 

-EE [(^"^" + c„r-"-i)(n6;C - (n + l)c^r-"-i)r 

n— k— — n 

00 n 

= E E l/nPN&nP(ro^"+^-rr+^)-(n + l)|c„p(ro 

n— k——n 
00 n 



n=0 k= — n 

The estimate (P?7)) yields 



00 n 



« ^E E "l/nl' (l^"^0"^' + Ic»^e-'"-^) ■ (2.17) 
n— k— — n 

(i) Suppose that ec = — e., = 1. With the notation in (j2.5p . we have 

Kn^ + n)(K^ - Kc)(k, - k™) - + 1)k^ + nKc)((n + 1)k„i + riK^)! w n^iS^ + p^"+^), 

and hence 

„2n ^ 2n 



It then follows from (PTZ)) that 

°° " <■ ^2n|yfc|2 
7i—0k— — n 

Let 



logp 

If n < Ns, then we have that 5 < p'"'. and hence 

^2 ^ „2n - ^2 _|_ „2n 

n=Ok=-n " n<Ns k=-n " 



TV, = i^. (2.20) 



k— — m \k— — m 
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for any m < Ns. By taking S to be n — 1,2,..., we see that if the fohowing holds 

n 

limsup(rero)"/^ ^ l/'l=oo, (2.21) 

n— J-oo , 

k— — n 

then there is a sequence {rik} such that 

lim S |„,| =00, (2.22) 

i.e., weak CALR occurs. 

Suppose that the source function / is supported inside the critical radius r» — yfr^T\, (and 
outside Te) and its Newtonian potential F cannot be extended harmonicaUy in < r*. Then 

/ n \ 1/" 

limsup ^ |/^-| >1/V^. (2.23) 

and consequently, (|2.2ip holds. This proves that if the source function / is supported inside the 
sphere of critical radius r*, then weak CALR occurs. 

If the source function / is supported outside the sphere of critical radius r.^ = ^rerp, then its 
Newtonian potential F can be extended harmonically in |a;| < r* + 277 for 77 > and 

E E 9t7^ ^ E E < cmi.,,,^^^^, < oo. (2.24) 

n— k—~n n— k— — n 

So is bounded regardless of (5 and CALR does not occur. 

Suppose that / is supported inside and [GCl] holds. Let {rij} be the sequence satisfying 



hm n,^M/„^j2 ^ oo. 



J^OO 

k— — nj 

If (5 = for some a, let j(q!) be the number in the sequence such that 
Then, we have 

* g2 _|_ „2n - ^ „2n 

n<JVi k=-n ' 7i<Ns k=-n " 

"3(a) 

2n,- 



as a ^ OO. So CALR takes place. This completes the proof of (i) 
To prove (ii) assume that Ec 7^ = 1- In this case, we have 

^2n J2n 



P 



and 



n=Ofc=-n ^ 



8 



The rest of proof of (ii) is the same as that for (i) . 
Suppose now that —Eg 7^ 1. If Cc = 1, then we have 



\bn\ 



P 



and 



dF 



dr 



Since the source function / is supported outside the radius r^, wc have 

OF 

jr- < C'||/||i2(R3), 

and Es is bounded independently of S. The case when Ec 7^ 1 can be treated similarly. 



A Gap property of dipoles 

In this appendix we show that the Newtonian potentials of dipole source functions satisfy the gap 
conditions [GCl] and [GC2]. We only prove [GCl] since the other can be proved similarly. 

Let / be a dipole in B,-, \ B^, i.e., /(x) = a • V(5y(x) for a vector a and y S Br, \ B^- Then 
its Newtonian potential is given by i^(x) = —a • VyG(x — y). It is well-known (see, for example, 
[13) ) that the fundamental solution G(x — y) admits the following expansion if |y| > |x|: 



00 n 



G(x-y) = -E E 



n— k——n 



2n+ 1 



n+l ■ 



So we have 



and hence 



00 n y ^ 



n=0 fc=- 



1 



2n+ 1 



-a - V 



|y|ri+l ^n\y ) 



We show that 



E -^00 as 00, 



(A.l) 
(A.2) 



and hence [GCl] holds. The following lemma is needed. 



Lemma A.l For any a and y on and for any positive integer n there is a homogeneous 
harmonic polynomial h of degree n such that 



and 



a-Vhif)^! 



max /i X < . 

|3c|=i n 



(A.3) 
(A.4) 
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Proof. After rotation if necessary, we may assume that y = (1, 0, 0). We introduce three homoge- 
neous harmonic polynomials of degree n: 

hi{x) := ^ [{xi + ix^T + {^1 - ^X2T\ , 
Zn 

/i2(x) := [(a;i + ^^2)" - {xi - ix2T] , 
Zm 

/i3(x) := [(a;i + ixs)" - {xi - ixs)"] . 
Zm 

Then one can easily see that 

V/ii(y) = (1,0,0), V/i2(y) = (0,1,0), V/i3(y) = (0,0,1). 

So if we define 

h = aihi + a2h2 + 03/13, 

then (|A.3() holds. 

Since ^ 

max |/i., (x)| < — for 7 = 1, 2, 3, 

|x|=i ' ■'^ ^' ~ n 

we obtain (|A.4|) using the Cauchy- Schwartz inequality. This completes the proof. □ 
Let a and y be two unit vectors, and let ft, be a homogeneous harmonic polynomial of degree 
n satisfying (jA.3|) and (|A.4p . Then h can be expressed as 

n 

/i(x) - «fc|xr>;f(S), 

k— — n 

where 

«fc = ^ / ft(X)y„^(X)d5. (A.5) 

Because of (|A.3p , we have 

i = a.vMy)< KI|a-v(|xr%f(S))|. 

A;— — n 

So there is say /cn, between —n and n such that 

KJ|a.V(|xri;f"(X))|>^-L_. (A.6) 
and (jA.Sp . it follows by using Jensi 
afe„P<^^JMx)ni;f"(X)pd^< 



On the other hand, from (jA.4[) and (jA.Sp . it follows by using Jensen's inequality that 

3 



Thus wc have 



for some C independent of n. 

Now one can see from (|A.ip that 



n|y|"+i 



for some C independent of n. Since |y| < r*, we obtain that 

E ^^M/'l' > ^^M/M' ^ ^ (^j ^ °« as n ^ 00, 

as desired. It is worth mentioning that the constants C in the above may be different at each 
occurrence, but are independent of n. 
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